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The energy spectrum of nuclear recoils in Weakly Interacting Massive Particle (WIMP) direct 
detection experiments depends on the underlying WIMP mass (strongly for light WIMPs, weakly 
for heavy WIMPs). We discuss how the accuracy with which the WIMP mass could be de- 
termined by a single direct detection experiment depends on the detector configuration and the 
WIMP properties. In particular we examine the effects of varying the underlying WIMP mass, 
the detector target nucleus, exposure, energy threshold and maximum energy, the local velocity 
distribution and the background event rate and spectrum. 
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1. Introduction 

The direct detection of WIMPs in the lab would not only directly confirm the existence of dark 
matter but would also allow us to probe the WIMP properties, in particular its mass. This would 
shed light on its nature and probe extensions of the standard model of particle physics. Furthermore 
definitive detection of the WIMP may well require consistent signals (i.e. with the same inferred 
WIMP properties) from direct detection, indirect detection and collider experiments. Here we 
give a brief overview of recent work [p on determining the WIMP mass from direct detection 
experiments (see also Refs. [§, [|, |]]). 

The differential event rate (number of events per unit energy, time and detector mass) has a 
roughly exponential energy dependence: [||] 
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where c\ and C2 are fitting parameters of order unity (which depend on the target mass number and 
energy threshold), (dR/dE)o is the event rate in the E — ► OkeV limit and F(E) is the form factor. 
The characteristic energy scale of the exponential, £r, depends on the WIMP mass, m x , and is 
given by 
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where is the target nuclei mass and v c is the local circular speed. For light WIMPs (m x <c wja) 
£r « raj, while for heavy WIMPs (m x S> tka) £r ~ const. In other words, for light WIMPs the 
energy spectrum is strongly dependent on the WIMP mass while for heavy WIMPs the dependence 
on the WIMP mass is far weaker. Consequently it should be easier to measure the mass of light 
(compared with the target nuclei) WIMPs than heavy WIMPs. 



2. Monte Carlo simulations 

We have used Monte Carlo simulations to examine how well a SuperCDMS like detector [j|] 
could determine the WIMP mass from the energies of observed WIMP nuclear recoil events. Our 
benchmark detector is composed of Ge, has a nuclear recoil energy threshold of E& = lOkeV and 
has no upper limit on the recoil energy. We assume that the detection efficiency is independent 
of energy, the energy resolution is perfect, the background is zero, the form factor has the Helm 
form and fix the WIMP-proton cross-section to be a p = 10~ 8 pb, a factor of a few below the current 
exclusion limits [g]. We assume the local WIMP speed distribution is Maxwellian and the local 
density is 0.3GeVcm -3 and consider (efficiency weighted) exposures of £ = 3 x 10 3 , 3 x 10 4 and 
3 x 10 5 kgday which correspond, roughly, to a detector with mass equal to that of the 3 proposed 
phases of SuperCDMS taking data for a year with a ~ 50% detection efficiency. These are, gen- 
erally, optimistic assumptions and will therefore give 'best case' results. We discuss the effects of 
dropping or varying some of these assumptions in Sec. ||, see also Ref. [[I]]. 

For each WIMP mass, m™, and detector configuration we calculate the probability distribution 
of the maximum likelihood estimator of the WIMP mass by simulating 10 4 experiments. We first 
calculate the expected number of events, A, from the input energy spectrum. The actual number of 
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Figure 1: The fractional deviation of the WIMP mass limits from the input mass, (m^ m — m™)/m™, for 
exposures c? = 3xl0 3 ,3xl0 4 and 3 x 10 5 kgday and input cross-section (7 p = 10~ 8 pb for the benchmark 
SuperCDMS like detector. The solid (dotted) lines are the 95% (68%) confidence limits. 



events for a given experiment, A^ eX pt» is drawn from a Poisson distribution with mean A. We Monte 
Carlo generate A'expt events from the input energy spectrum, from which the maximum likelihood 
mass and cross-section for that experiment are calculated. Finally we find the (two-sided) 68% and 
95% confidence limits on the WIMP mass from the maximum likelihood masses. 



3. Results and discussion 

The accuracy with which the WIMP mass could be measured by the benchmark SuperCDMS ^ 
like Ge detector described above is shown in Fig. [l[ With exposures of <§ = 3 x 10 4 and 3 x 
10 5 kgday it would be possible to measure the mass of a light, m x ~ <^(50GeV), WIMP with an 
accuracy of roughly 25% and 10% respectively. For heavy WIMPs (m x 3> lOOGeV) even with a 
large exposure it will only be possible to place a lower limit on the mass. For very light WIMPs, 
m x < ^(20GeV), the number of events above the detector energy threshold would be too small to 
allow the mass to be measured accurately. 

The number of events detected is directly proportional to both the exposure and the cross- 
section, therefore these quantities have the greatest bearing on the accuracy of the WIMP mass 
determination. 

The energy threshold, and the maximum energy, E max , above which recoils are not de- 
tected/analysed also affect the accuracy with which the WIMP mass can be determined. Increasing 
.Eth ( or decreasing £m ax ) not only reduces the number of events detected, but also reduces the range 
of recoil energies and the accuracy with which the characteristic energy of the energy spectrum, 
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£r, and hence the WIMP mass, can be measured. For light WIMPs the small £r means that the 
expected number of events decreases rapidly as the energy threshold is increased, while for heavy 
WIMPs the large £r, and flatter energy spectrum, means that the smaller range of recoil energies 
reduces the accuracy with which £r can be measured. Reducing the maximum energy only has a 
significant effect for heavy WIMPs. 

The WIMP and target mass dependence of Z?r suggests that heavy targets will be able to 
measure the mass of a heavy WIMP more accurately, however the rapid decrease of the nuclear 
form factor with increasing momentum transfer which occurs for heavy nuclei means that this is in 
fact not the case (see also Ref. [Q]). 

If the WIMP distribution on the ultra-local scales probed by direct detection experiments is 
smooth, then the ±20kms _1 uncertainty in the local circular speed [[7|] leads to a ~ 10% systematic 
error in the determination of m x . Changes in the detailed shape of the local velocity distribution 
lead to relatively small changes in the shape of the differential event rate and hence a relatively 
small, &(5%), systematic uncertainty in the WIMP mass. If the ultra-local WIMP distribution 
consists of a finite number of streams, then the energy spectrum will consists of a number of steps. 
The positions of the steps will depend on the (unknown) stream velocities, as well as the target 
nuclei and WIMP masses. With multiple targets it would in principle be possible to constrain the 
WIMP mass without making any assumptions about the WIMP velocity distribution 

Future experiments aim to have negligible backgrounds, however if the background rate is not 
negligible compared with the WIMP event rate it will be difficult to disentangle a WIMP signal 
(and the WIMP mass) from the background if the background spectrum has a similar shape to 
the WIMP spectrum (i.e. exponential background, or flat background and a heavy WIMP). The 
uncertainties from backgrounds could be mitigated by using multiple targets and/or using multiple 
scatter events to measure/constrain the background spectrum. 
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